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Abstract 

An extensive program for the calculation of galactic cosmic -ray propagation has been de- 
veloped. Primary and secondary nucleons, primary and secondary electrons, secondary 
positrons and antiprotons are included. Fragmentation and energy losses are computed 
using realistic distributions for the interstellar gas and radiation fields. 

Models with diffusion and convection only do not account naturally for the observed 
energy dependence of B/C, while models with reacceleration reproduce this easily. The 
height of the halo propagation region is determined, using recent 10 Be/ 9 Be measure- 
ments, as greater than 4 kpc. The radial distribution of cosmic-ray sources required is 
broader than current estimates of the SNR distribution for all halo sizes. Our results in- 
clude an estimate of cosmic-ray antiproton and positron spectra, and the Galactic diffuse 
7-ray emission (see accompanying paper: Moskalenko 1998b). 
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Introduction. 



We are constructing a model which aims to reproduce self-consistently observational 
data of many kinds related to cosmic-ray (CR) origin and propagation: direct measure- 
ments of nuclei, electrons, positrons, antiprotons, gamma rays, and synchrotron radia- 
tion. These data provide many independent constraints on any model and our approach 
is able to take advantage of this since it must be consistent with all types of observation. 
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Here we present our results on the evaluation of diffusion/convection and reacceler- 
ation models based on the B/C and 10 Be/ 9 Be ratios, and set limits on the halo size. A 
re-evaluation of the halo size is desirable since new 10 Be/ 9 Be data are now available 
from Ulysses (Connell 1998) with better statistics than previously. Our preliminary re- 
sults were presented in Strong (1997a,b) and full results for protons, Helium, positrons, 
and electrons in Moskalenko (1998a). Some illustrative results for gamma-rays and syn- 
chrotron radiation are given in Strong (1997a) and Moskalenko (1998b) and all details 



are given in Strong (1998^1 
The model description. 

The models are three dimensional with cylindrical symmetry in the Galaxy, and the ba- 
sic coordinates are (R, z, p), where R is Galactocentric radius, z is the distance from the 
Galactic plane, and p is the particle momentum. The propagation equations are solved 
numerically on a grid by the method described in Strong (1998). R & is taken as 8.5 kpc. 
The propagation region is bounded by R = R h , z = Zh beyond which free escape is 
assumed. We take R h = 30 kpc. The range z h = 1 — 20 kpc is considered. For a given Zh 
the diffusion coefficient as a function of momentum is determined by B/C for the case 
of no reacceleration; if reacceleration is assumed then the reacceleration strength (re- 



'For more details see http://www.gamma.mpe- garching.mpg.de/^aws/aws.html 



05 



0.40 



0.30 



0.20 - 



0.10 - 



R = 8.5 

z — 0.0 kpc 



model 

15-003000 
\\ 15-003010 
\ 15-003020 



0.00 




03 



0.40 


i 


i i ; 

R = 8.5 

z — 0.0 kpc 


0.30 




n \ 

\ model 

\\ \l \ 15-005501 
\\ \ 1 5 — 005504 
\\| \ 15-005505 
1 \ 15-005507 


0.20 




' \\ \ 
•\'\ 




■■// 

/,// 


0.10 




"\\ 


0.00 


i 


I I • 



10' 10^ 10 u 10^ 10 u 10" 

kinetic energy. MeV/nucleon 



10' 10* 10' 10^ 10 J 10 L 

kinetic energy, MeV/nvccleon 



Fig. 1. Left panel: B/C ratio for diffusion/convection models without break in diffusion 
coefficient (Strong 1998 and references therein), for = 3 kpc, dV/dz = (solid line), 5 (dotted 
line), and 10 km s _1 kpc -1 (dashed line); solid line: interstellar ratio, shaded area: modulated 
to 300 - 500 MV; data: HEAO-3, Voyager, Ulysses. Right panel: B/C ratio for diffusive 
reacceleration models (Strong 1998) with zh = 5 kpc, va = (dotted), 15 (dashed), 20 (thin 
solid), 30 km s _1 (thick solid). In each case the interstellar ratio and the ratio modulated to 500 
MV is shown. 

lated to the Alfven speed) is constrained by the energy-dependence of B/C. The spatial 
diffusion coefficient for the case of no reacceleration is taken as D xx = (3D (p/ p ) Sl be- 
low rigidity p , f3D (p/ p ) 5 ' 2 above rigidity p . Since the introduction of a sharp break in 
D xx is an extremely contrived procedure which is adopted just to fit B /C at all energies, 
we also consider the case S± = S 2 (no break). The convection velocity V = V(z) is as- 
sumed to increase linearly with distance from the plane. For the case with reacceleration, 
the spatial diffusion coefficient is D xx = (3D (p/p ) 5 with 5 = | for all rigidities, and 
the momentum-space diffusion coefficient D pp is related to D xx (using Berezinskii 1990 
and Seo 1994). The injection spectrum of nucleons is assumed to be a power law in mo- 
mentum. The interstellar hydrogen distribution uses HI and CO surveys and information 
on the ionized component; the Helium fraction of the gas is taken as 0.11 by number. 
Energy losses for nucleons by ionization and Coulomb interactions are included. The 
distribution of CR sources is chosen to reproduce the CR distribution determined by 
analysis of EGRET 7-ray data (Strong 1996b). The secondary nucleon and secondary 
e ± source functions are computed from the propagated primary distribution and the gas 
distribution. 



Illustrative results. 

We consider the cases of diffusion+convection and diffusion+reacceleration, since these 
are the minimum combinations which can reproduce the key observations. Our basic 
conclusion is that the reacceleration models are more satisfactory in meeting the con- 
straints provided by the data, reproducing the B/C energy dependence without ad hoc 
variations in the diffusion coefficient; further it is not possible to find any simple version 
of the diffusion/convection model which reproduces B/C satisfactorily. 

Figure la shows the diffusion+convection model without break, 5\ = 62', for each 
dV/dz, the remaining parameters D , Si and p are adjusted to fit the data as well as 
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Fig. 2. (a) Predicted 10 Be/ 9 Be ratio as function of 2^ for dV/dz = 0, 5, 10 km s _1 kpc -1 , the 
Ulysses experimental limits are shown as horizontal dashed lines. The shaded regions show the 
parameter ranges allowed by the data, (b) 10 Be/ 9 Be ratio for diffusive reacceleration models as 
function of at 525 MeV/nucleon corresponding to the mean interstellar value for the Ulysses 
data (Connell 1998). 

possible. It is clear that a good fit is not possible; the basic effect of convection is to re- 
duce the variation of B/C with energy, and although this improves the fit at low energies 
the characteristic peaked shape of the measured B/C cannot be reproduced. If we allow 
81 7^ 5 2 it can clearly be fitted, but the break has to be large and the procedure is ad hoc. 

Figure lb illustrates a diffusive reacceleration model and shows the effect on B/C of 
varying va (— 30 km s -1 ) for z h = 5 kpc. This shows how the initial form becomes 
modified to produce the characteristic peaked shape. Reacceleration models thus lead 
naturally to the observed peaked form of B/C, as pointed out by previous authors (e.g., 
Letaw 1993, Seo 1994, Heinbach 1995); a value va ~ 20 km s -1 seems satisfactory. 

Figure 2a summarizes the limits on z h and dV/dz for diffusion/convection, using 
the 10 Be/ 9 Be ratio at the interstellar energy of 525 MeV/nucleon appropriate to the 
Ulysses data (Connell 1998). We conclude that in the absence of convection 4 kpc < 
Zh < 12 kpc, and if convection is allowed the lower limit remains but no upper limit 
can be set. In the case dV/dz < 7 km s -1 kpc -1 , this figure places upper limits on the 
convection parameter for each halo size. These limits are rather strict, and a finite wind 
velocity is only allowed in any case for z h > 4 kpc. Figure 2b shows 10 Be/ 9 Be for 
the reacceleration models as a function of z h at 525 MeV/nucleon corresponding to the 
Ulysses measurement and we again find that 4 kpc < Zh < 12 kpc. 

Figure 3 (left panel) shows the effect of halo size on the radial distribution of 3 
GeV CR protons, for the reacceleration model. For comparison we show the CR distri- 
bution deduced by model-fitting to EGRET gamma-ray data (> 100 MeV) from Strong 
(1996b), which is dominated by the 7r°-decay component; the analysis by Hunter (1997), 
based on a different approach, gives a similar result. The predicted CR distribution using 
the SNR source function is too steep even for large halo sizes; in fact the halo size has a 
relatively small effect on the distribution. Other related distributions such as pulsars have 
an even steeper falloff. Based on these results we have to conclude, in the context of the 
present models, that the distribution of sources is not that expected from the (highly un- 
certain) distribution of SNR. In view of the difficulty of deriving the SNR distribution 
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Fig. 3. Left panel: radial distribution of 3 GeV protons at z = 0, for diffusive reacceleration 
model with halo sizes Zh = 1, 3, 5, 10, 15, and 20 kpc (solid curves). Dashed line: the source 
distribution is that for SNR given by Case (1996), histogram: the CR distribution deduced 
from EGRET >100 MeV gamma rays (Strong 1996b). Right panel: radial distribution of 3 
GeV protons at z = for the source distribution actually adopted (dashed line), for diffusive 
reacceleration model with various halo sizes z^ = 1, 3, 5, 10, 15, and 20 kpc (solid curves). 

this is perhaps not a serious shortcoming; if SNR are indeed CR sources then it is pos- 
sible that the gamma-ray analysis gives the best estimate of their Galactic distribution. 
Therefore, we have chosen a CR source distribution to fit the 7-ray data after propaga- 
tion (Figure 3, right panel). The possibility of anisotropic diffusion (preferentially in the 
radial direction) has not yet been addressed in our models. 

The positron fraction computed is in good agreement with the measured one be- 
tween 1 and 10 GeV, where the data are rather precise. Our positron predictions from 
Moskalenko (1998a) have been compared with more recent absolute measurements in 
Barwick (1998) and the agreement is good; for the positrons this new comparison has the 
advantage of being independent of the electron spectrum (see also Moskalenko 1998b). 
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